The effective force NL3 revisited 
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Covariant density functional theory based on the relativistic mean field (RMF) Lagrangian with 
the parameter set NL3 has been used in the last ten years with great success. Now we propose a 
modification of this parameter set, which improves the description of the ground state properties of 
many nuclei and simultaneously provides an excellent description of excited states with collective 
character in spherical as well as in deformed nuclei. 
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Density functional theory is a universal and power- 
ful tool for describing properties of finite nuclei all over 
the periodic table. In the non-relativistic framework the 
most successful density functionals are the ones based on 
density dependent forces, such as the Skyrme [l| or the 
Gogny 2] functional. Relativistic mean field (RMF) the- 
ory was first introduced as a fully fledged quantum field 
theory by Walecka [!, Q|. However, it turned out very 
soon [f|, that for a quantitative description of nuclear 
surface properties an additional density dependence is 
necessary. Nowadays RMF theory modified in this form 
is considered as a covariant form of density functional 
theory. Over the years it has gained considerable interest, 
in particular, for the description of nuclei at and far from 
stability [f| @, H, G3 • Compared with non-relativistic 
density functionals covariant density functional theory 
has certain advantages. They are characterized by a new 
saturation mechanism obtained by a delicate balance be- 
tween a strongly attractive scalar field and a strongly re- 
pulsive vector field. Moreover, the very large spin-orbit 
splitting, observed in finite nuclei, is a relativistic effect. 
Therefore, its treatment in relativistic models arises in a 
natural way without any additional adjustable parame- 
ters. In addition, time-odd mean fields which are impor- 
tant in systems with broken time reversal symmetry are 
uniquely defined in RMF theory because of the Lorentz 
covariance of the underlying Lagrangian (Tlj . 

Pairing properties are essential for a description of nu- 
clei with open shells. They have been included first in the 
constant gap approximation by occupation numbers of 
BCS-type [12j. Since this procedure requires the knowl- 
edge of the experimental pairing gaps, it cannot be ap- 
plied in unexplored regions of the nuclear chart, where 
the binding energies are not known. In addition, it is 
noted that the BCS approximation breaks down in nu- 
clei far from the valley of stability, where the coupling 
to the continuum is essential [l3| . Therefore the con- 
stant gap approximation has been replaced by relativistic 
Hartree-Bogoliubov (RHB) theory [lj] which includes a 
finite range particle-particle interaction of Gogny form. 
The details of this theory have been discussed in several 
review articles @, 0, [1] and in the references given there. 

In any case the adopted functionals are considered uni- 



versal in the sense that they can be used for nuclei all 
over the periodic table, where mean field theory is ap- 
plicable. It is therefore very desirable to find a unique 
parameterization for the Lagrangian of the model, which 
is able to describe as many experimental data as pos- 
sible. In other words, we search for an effective force 
that is able to describe properties of nuclei from light 
to very heavy, from the proton to the neutron drip line. 
Moreover, a powerful density functional should not only 
describe the ground state properties of finite nuclei but 
also, at the same time, collective excited states within 
time-dependent density functional theory. This is crucial 
because otherwise, one would end up with one effective 
force for the description of the masses, another one for 
exotic systems, a third one for giant resonances an so on. 
Definitely, such specific parameterizations do not serve 
the purposes of universal density functional theory. Their 
inability originates from the fact that they are tailored 
to describe either specific observables or narrow regions 
of the periodic chart. 

There are several types of relativistic density function- 
als. Conventional RMF theory is based on the Walecka 
model, where the nucleons interact by the exchange of 
phenomenological mesons and, up to now, in all the suc- 
cessful relativistic models two additional assumptions are 
essential: (i) in the mean field approximation only the 
nucleonic fields are quantized and the nucleons move in- 
dependently in classical mesons fields, depending in a 
self-consistent way on the nuclear densities and currents, 
and (ii) in the no-sea approximation vacuum polarization 
and the contributions of the negative energy solutions are 
not explicitly taken into account. 

Of course, the pure Walecka model is unable to de- 
scribe nuclear matter and finite nuclei in a quantitative 
way. The medium dependence induced by many-body 
correlations is not reproduced properly. This leads to an 
effective density dependence, which can be taken into ac- 
count either by non-linear self-interactions between the 
mesons or by a density dependence of the coupling con- 
stants. Models with density dependent coupling con- 
stants are tailored more in the spirit of density functional 
theory. Indeed, they have turned out to be very suc- 
cessful [ill EH G3- However, they are technically more 
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TABLE I: The experimental values of the energies of the Giant Monopole (GMR) and Giant Dipole (GDR) Resonances in 
208 Pb are compared with the predictions of various non-linear forces. 





Exp. 


NL3 


NL3* 


NL-SH 


NL-1 


NL-Z 


NL-RA1 


GMR 


14.1 ±0.3 


14.04 


13.90 


16.48 


12.44 


11.72 


15.26 


GDR 


13.3 ±0.1 


12.95 


12.95 


12.98 


12.97 


12.53 


12.87 



complicated, in particular, for RPA calculations. On 
the other hand, non-linear models have some justifica- 
tion. In fact, the effects of vacuum polarization which 
can contribute up to 30 % to the resulting binding en- 
ergies are not taken into account explicitly in the frame- 
work of the no-sea approximation, but only globally by 
adjusting the parameters of the model. A proper renor- 
malization procedure in nuclear matter (l8j and in finite 
nuclei [HI, [13, H3 leads to counterterms in the form of 
non-linear meson coupling. From this point of view non- 
linear meson coupling theories seem at a first glance to be 
more appropriate for a phcnomenological description of 
vacuum polarization as compared to those using density 
dependent coupling constants. Of course, by adjustment 
of the final Lagrangian to the same experimental data, 
both methods to introduce a density dependence lead to 
very similar results. However the details are different: 
non-linear meson couplings lead in the case of nuclear 
matter to density dependent meson masses, i.e. to den- 
sity dependent range parameters of the force and no rear- 
rangement terms in the Dirac equation, whereas density 
dependent coupling constants lead to density dependent 
strenght parameters of the effective force and additional 
rearrangement terms. So far such differences have not 
been investigated in detail. 

The parameter set NL3 (22| with non-linear meson cou- 
plings represents one of most successful non-linear RMF 
forces. It was proposed ten years ago. In the meantime, 
new experimental data on nuclear masses have appeared. 
Moreover, new and more reliable information about the 
neutron skin became available. On the other hand, it was 
found that NL3 encounters some difficulties in describing 
light Hg and Pb isotopes (23[ and certainly, there is al- 
ways a need for better predictions of the masses which, of 
course, reflect also correct nuclear sizes. For this reason 
we decided to improve the parameter set NL3 by perform- 
ing a new global fit of ground state properties of spherical 
nuclei and infinite nuclear matter. New parameterization 
obtained in this fit will be called NL3*. 

In order to stay as close as possible to the very suc- 
cessful parameter set NL3, we kept in this fit the func- 
tional form of the Lagrangian unchanged. In particular 
the density dependence of the both parameter sets NL3 
and NL3* is determined by the non-linear meson cou- 
plings in the isoscalar channel a. We will discuss the 
adequacy of this ansatz at the end the this letter. 

The starting point of Covariant Density Functional 



Theory (CDFT) is a standard Lagrangian density (24| 

C = 4> (<y(id - g u uj - g p pf - eA) - m - g a a) ip 

+\{daf -\rnlo 2 -\^ V W» + \mW (1) 

which contains nucleons described by the Dirac spinors ip 
with the mass m and several relativistic fields character- 
ized by the quantum numbers of spin, parity, and isospin. 
These arc effective fields mediated by mesons, with no di- 
rect connection to mesons and resonances existing in free 
space. They only carry their quantum numbers and char- 
acterize the properties of the possible relativistic fields 
entering the effective Dirac e quat ion, which corresponds 
to the Kohn-Sham equation [25j | in the non-relativistic 
case. It is only for simplicity that we use the conventional 
names a {r = 0+,T = 0), u (I* = 1~,T = 0), and p 
(J* = l~,T = 1). In addition, we have the electromag- 
netic field A. In principle there should also be an effective 
scalar meson with I" = + , and T = 0, the 5-meson. It 
plays an important role in G-matrix calculations [26| . 
However, present phenomenological parameterizations of 
the relativistic energy density functional 0, [l?], HH are 
very successful without the inclusion of this meson. On 
the basis of present experimental data in finite nuclei it is 
very difficult to distinguish between effective scalar and 
vector mesons in the isospin T = 1 channel, and, there- 
fore, all the isospin dependence is carried by the effective 
p-meson. 

The Lagrangian |T]) contains as parameters the me- 
son masses m a , m w , and m p and the coupling constants 
5<t, fl'wj an d SL>- This model has first been introduced by 
Walecka jH, 0| . It soon became clear that surface prop- 
erties of finite nuclei, in particular, the incompressibility, 
cannot be described properly by this model. Therefore, 
Boguta and Bodmer [j| introduced a non-linear meson 
coupling 

U(a) = \mla 2 + ^g 2 a z + - A g^ A . (2) 

which brings in an additional density dependence. By 
fitting experimental data of some carefully chosen spher- 
ical nuclei, very successful parameter sets have been pro- 
posed. For example, the set NL1 [l^] is an excellent force 
for describing ground state and rotational properties of 
finite nuclei along the valley of stability. However, it pro- 
duces a very large asymmetry parameter 04. Therefore, 
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TABLE II: The total binding energies B.E., charge radii r c , 
and the differences between the radii of neutron and proton 
density distributions r np = (r n — r p ), used to adjust the new 
parameter set NL3*. The calculated values are compared with 
experimental data (values in parentheses). 



B.E.(MeV) 



r c (fm) 



r„-r p (fm) 



16Q 

40 Ca 
48 Ca 
72 Ni 
90 Zr 
116 Sn 
124 Sn 
132 Sn 

204p b 
208p b 
214p b 
210p o 



-128.112( -127.619) 
-341.578( -342.052) 
-413.615( -415.990) 
-612.168( -613.152) 
-782.368( -783.891) 
-986.512( -988.680) 



2.735(2.730) -0.15 

3.470(3.485) -0.14 
3.470(3.484) 0.14 
3.892 0.26 
4.263(4.272) 0.04 (0.07) 
4.604(4.626) 0.10 (0.12) 



-1048.324(-1049.962) 4.655(4.674) 0.22 (0.19) 
-1101.550(-1102.850) 4.636 
-1608.100(-1607.505) 
-1638.230(-1636.430) 



0.29 

5.432(5.486) 0.19 
5.508(5.505) 0.23 (0.20) 



-1660.119(-1663.290) 5.566(5.562) 0.26 
-1649.024(-1645.210) 5.544 0.20 



later, the parameter set NL3 |22| was introduced to im- 
prove the asymmetry parameter using the same ansatz 
without increasing the number of phenomenological pa- 
rameters. This set is able to describe many ground state 
properties of finite nuclei all over the periodic table. Re- 
cently, the set NL3 has also been applied very successfully 
for the investigation of excited states such as collective ro- 
tations [13, [H, Hi| and giant resonances [3(| HH H3, [33[ . 
In particular, the use of the effective force NL3 in the 
p/i-channel together with the Brink-Booker part of the 
Gogny force D1S in the pp-channel [14j . within the frame- 
work of RHB theory, led to a very successful tool for the 
description of many properties of ground states as well 
as of excited states with collective character. The NL3 
force was for a long time the only parameterization of the 
Lagrangian of the RMF theory that was able to provide 
simultaneously a consistent description of both excitation 
energies of giant monopole and giant dipole resonances. 

In Table HI the experimental values of the energies 
of the Giant Monopole (GMR) and Giant Dipole Res- 
onances (GDR) in 208 Pb are compared with the predic- 
tions of various non- linear forces. It is clearly seen that all 
forces give reasonable values for the GDR energy, which 
are close to experiment. This is not the case, however, 
for the giant monopole resonance. There, all forces pre- 
dict values which, either overestimate (NL-SH 1341 and 
NL-RA1 [H) or underestimate (NL1 Qj, NLZ 0) the 
experimental value. Only NL3 and NL3* give results 
which are in excellent agreement with experiment. This, 
of course, is connected with the value of the incompress- 



TABLE III: Parameters of the effective interaction NL3* in 
the RMF theory together with the nuclear matter properties 
obtained with this effective force. The values of nuclear mat- 
ter properties obtained with NL3 are shown in parentheses. 

Parameters of NL3* 



M 


= 939 (MeV) 






rn a 


= 502.5742 (MeV) 


9<? 


= 10.0944 




= 782.600 (MeV) 


gu, 


= 12.8065 


m p 


= 763.000 (MeV) 


9p 


= 4.5748 


<n 


= -10.8093 (fm" 1 ) 






93 


= -30.1486 







Nuclear matter properties 



po = 0.150 
(E/A) x = 16.31 
K = 258.28 
J = 38.6 
m*/m = 0.594 



(0.148) fm" 3 
(16.30) MeV 
(271.76) MeV 
(37.4) MeV 
(0.60) 



ibility of nuclear matter predicted by each of these forces. 

The starting point for the fit of the new parameter set 
were the NL3 values, this time, however, the mass m u of 
the w-meson was kept fixed to to w = 782.6 MeV. This 
value differs slightly from that obtained in the NL3 fit 
in Ref. [22j and is the same as the one proposed for the 
free nucleon-nucleon force by Machleidt in Ref. [37j • As 
shown in Table |TT] we used twelve spherical nuclei for the 
present fit: 16 0, 40 Ca, 48 Ca, 72 Ni, 90 Zr, 116 Sn, 124 Sn, 
204 Pb, 208 Pb, 214 Pb, and 210 Po. Compared with the in- 
put for the NL3 fit in Ref. [13| we added the nuclei 204 Pb 
and 210 Po and replaced 58 Ni by 72 Ni in order to include 
more data on neutron rich isotopes. The experimental 
binding energies were taken from Ref. [38| , and the charge 
radii from a recent compilation in Ref. [3^ |. The most 
recent experimental information on the neutron skins of 
the nuclei 90 Zr 0o|, 116 Sn, 124 Sn and 208 Pb [MIH was 
also included. The values in parentheses correspond to 
the relative error bars used in the fitting procedure. The 
binding energies of finite nuclei and the charge radii are 
taken within an uncertainty of 0.1% and 0.2%, respec- 
tively. The error used for the neutron skin is 5%. 

For the open shell nuclei pairing correlations are 
treated in the BCS approximation with empirical pair- 
ing gaps (five-point formula). Nuclear matter informa- 
tion was also considered in the fit. The "empirical" in- 
put is: E/A = -16 MeV (5%), p = 0.153 fm" 3 (10%), 
A'o = 250 MeV (10%), and J = 33 MeV (10%). The 
values in parentheses correspond to the error bars used 
in the fitting procedure. The resulting values of the La- 
grangian parameterization NL3* are given in Table IIII1 
together with the predictions of the force for the nuclear 
matter properties at the saturation point. The values 
in parentheses are the predictions obtained with the ef- 
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FIG. 1: Absolute deviations of the binding energies calculated 
with the parameter sets NL3 (filled circles) and NL3* (open 
circles) from the experimental values of Ref. [38l ]. 



fective interaction NL3. The resulting energy at satura- 
tion and asymmetry energy are rather close to the former 
values, only incompressibility is about 12 MeV smaller. 
Nonetheless both values are still sowewhat to high as 
compared with other paremeter sets with a density de- 
pendence in the iso- vector channel [43|, [44[ ■ 

In Fig. [1] the binding energies of more than 180 even- 
even nuclei are compared with experiment and the pre- 
dictions of the NL3 forces. All calculations have been 
performed within the RHB model with the Gogny force 
D1S [451 ] in the pairing channel. 



V pp (l,2) 



i=l,2 



(Wi + B l P a -H l P T -M l P IJ P T ) 



(3) 

Hi, and Mi (i = 1,2) of 
. It should be emphasized 



The parameters Wi, B 
this force are given in Ref. [4 
that, since a finite-range pairing interaction is used, the 
results do not depend on unphysical parameters like, the 
momentum cut-off in the pairing channel. 

The results of NL3* are shown as open circles while 
those obtained with NL3 are marked by filled circles. For 
light nuclei, both forces give similar predictions, however, 
as the mass number increases NL3* results are clearly 
closer to the zero MeV line. 

As discussed above the parameter set NL3 had difficul- 
ties to reproduce the proper ground state deformations 
in light Hg and Pb nuclei [23|]. This is no longer the 
case with the parameter set NL3*. To investigate this, 
we have carried out constrained axially deformed RHB 
calculations of several even-A Pb isotopes with masses 
between 182 < A < 192 in an external quadrupole field 
and we display in Fig. [3 the corresponding energy sur- 
faces as a function of the quadrupole deformation. It 
is seen that in all cases the Pb isotopes turn out to be 
spherical. This is also the case for all other Pb isotopes, 
which are not shown in the figure. It is also seen that 
the Pb isotopes manifest the interesting effect of shape 



coexistence. The energies which correspond to the oblate 
and prolate shape solutions are very close to the spherical 
ones but definitely lay higher in energy. This indicates 
a clear improvement as compared to the parameter set 
NL3, where some light Pb isotopes showed a deformed 
shape 23]. A more quantitative analysis goes beyond 
the mean field limit and requires, for instance, GCM- 
calculations [47|, |48(. It is essential, however, that the 
mean field solution, which is the starting point for such 
investigations gives the correct behavior. 

In the following we investigate dynamical processes 
such as collective vibrations with the same parameter 
set NL3*. For that purpose we study the time-dependent 
RMF or RHB equations in the small amplitude limit [Zi| , 
i.e. we solve the relativistic RPA or QRPA equations. In 
Fig. [3] we display results for the monopole and isovector 
dipole response for the nucleus 208 Pb. For the multipole 
operator the response function R{E) is defined 



R(E) = ]T B{J t - 7 ) 



r/27r 



(£-£ 4 ) 2 + r 2 /4' 



(4) 



where Y is the width of the Lorcntzian distribution, and 



B{J% - Of) 



1 



2Ji + 1 



K^IIQaIIO/)! 5 



(•5) 



In the examples considered here the continuous strength 
distributions are obtained by folding the discrete spec- 
trum of R(Q)RPA states with the Lorentzian (see Eq. 
(3D) with constant width Y = 0.5 MeV (32|. The cal- 
culated peak energies of the ISGMR resonance at 13.9 
MeV and of the IVGDR resonance at 12.95 MeV should 
be compared with the experimental excitation energies: 
E = 14.1 ± 0.3 MeV [5(| for the monopole resonance, 
and E = 13.3±0.1 MeV [5l[ for the dipole resonance, 
respectively. Clearly, the agreement with experiment is 
very good. The NL3 force predicts the same value for 
the IVGDR and a slightly larger value for the ISGMR. 
This is to be expected as the predicted value for the nu- 
clear incompressibility with NL3 is larger than the one 
obtained with NL3*. However, both values are within 
the experimental error bars and one could say that RMF 
theory with non-linear meson couplings leads to values 
of the incompressibility around 260 ±10 MeV. 

Recently a new computer code has been developed for 
the solution of the relativistic QRPA equations in axially 
deformed nuclei [H, US ■ We used this code for the study 
of giant resonances in deformed nuclei. Here, we present 
as an example calculations in the prolate deformed nu- 
cleus 100 Mo. In Fig. 0] we show the total isovector dipole 
cross section as function of the GDR energy. The param- 
eter set NL3* is very effective in reproducing these exper- 
imental data. The full line corresponds to the fully self- 
consistent deformed relativistic QRPA calculations while 
the dotted line are the experimental data 0, [H| . The 
estimated centroid energy for the GDR differs from the 
experimental value by less than 0.2 MeV. It is noted that 
NL3 also predicts excellent results, however, our analysis 
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p, (3, FIG. 4: The total isovector electric dipole cross section (in 

mb) of the deformed nucleus 100 Mo as a function of the ex- 
citation energy, calculated with the effective force NL3*, is 
FIG. 2: Excitations energies of even-Pb isotopes as a function compared with experiment [H, H . 
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FIG. 3: The isoscalar monopole (a), and the isovector dipole 
(b) strength distributions in 208 Pb calculated with the effec- 
tive interaction NL3*. The experimental excitation energies 
are 14.1 ± 0.3 MeV [5(|] for the monopole resonance, and 
13.3 ±0.1 MeV [HH for the dipole resonance, respectively. 



shows that the results with the newly developed NL3* 
are slightly better. This can be traced back to the im- 
provement in the density dependence of the non-linear 
sigma channel of the new force. 

The response of the nuclei to an external force pro- 
vides an important test of our understanding of their 
microscopic structure. The Coriolis force is one of such 
forces. It is of particular interest because it breaks time 
reversal symmetry in the intrinsic frame and leads there- 
fore to currents and to nuclear magnetic fields having its 
origin in the spatial contributions of the vector mesons 
[56(. The impact of this field can be studied in rotat- 



ing nuclei. Such studies were previously performed with 
various non-linear RMF parameterizations in the differ- 
ent parts of the nuclear chart, see Ref. [7J and references 
therein. In order to test the parameter set NL3* we se- 
lected several representative examples. This choice was 
guided by (i) the necessity to test the accuracy of the 
description of different types of rotational bands in dif- 
ferent mass regions, and (ii) the restriction to consider 
only bands which are very weakly affected by the pair- 
ing correlations. The later constraint is related to the 
fact that the pairing modifies considerably the rotational 
properties and some inadequacies in the pairing channel 
may prevent meaningful conclusions about the reliability 
of the parameterization. For a discussion of this point 
see Ref. [28 1. Thus, a highly-deformed band in the nu- 
cleus 58 Cu 57] . a non-terminating ground state band in 
74 Kr [58|, a smoothly- terminating band in 109 Sb [59] . 
and an yrast superdeformed (SD) band in 143 Eu [6(| were 
selected for this comparison. The pairing correlations are 
weak in these bands, and thus all calculations were per- 
formed within the framework of cranked relativistic mean 
field (CRMF) theory [6l| which neglects pairing. 

The results of these calculations for rotational and de- 
formation properties of the bands mentioned before are 
compared with experiment in Fig. [5] The kinematic 
(J*- 1 -*) and dynamic (J*- 2 -*) moments of inertia (Fig. 5a) 
of the highly-deformed band in 58 Cu as well as its transi- 
tion quadrupole moment Qt (Fig. 5e) are well described 
by the CRMF calculations. The fact that the jW mo- 
ment is considerably larger than the J^ 2 ' moment and 
that both moments smoothly decrease with rotational 
frequency are clear fingerprints that the pairing is not 
important at high spin [7| . While this condition is clearly 
satisfied in 58 Cu (see Fig. 5a), these moments are closer 
to each other in the case of the SD band in 143 Eu (Fig. 
5b) suggesting that the pairing is relatively more impor- 
tant in the later band. As a result, we observe somewhat 
larger discrepancies between experiment and theory in 
the case of 143 Eu (Fig. 5b). If the calculations would 
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Rotational frequency £2 X [MeV] Spin I Spin I 

FIG. 5: Panels (a-b):. Kinematic (solid circles, solid lines) and dynamic (open circles, dashed lines) moments of inertia in 
the highly-deformed band in 58 Cu and the yrast superdeformed band of 143 Eu. Experimental data are shown by circles, while 
the results of calculations by lines. Panels (c-d): Excitation energies of the non-terminating band in 74 Kr and the smoothly 
terminating band in 109 Sb relative to a rigid rotor reference Erld{I) = AI(I + 1), where A is a moment of inertia parameter. 
Panels (e-h): Experimental and calculated transition quadrupole moments Q t . Only the data for 143 Eu include the errors due 
to the uncertainties in the nuclear and electronic stopping powers. Other data exclude these uncertainties which are on the 
level of 10-15% of absolute value of Qt. The box in panel (e) indicates the upper and lower limits of Qt and the spin range 
where they have been measured. The panel (h) itself acts as such a box in the case of the 143 Eu band. Experimental data are 
taken from Refs. ( 58 Cu), ||| ( 74 Kr), [HEl] ( 109 Sb), and [6(j ( 143 Eu). 



include weak pairing, this would lead to the decrease (in- 
crease) of calculated J' 1 ' ( J^ 2 ') moments of inertia bring- 
ing them closer to experimental data. It is also necessary 
to mention that the SD band in 143 Eu is not linked to 
the low-spin spectra. 

Based on the effective alignment approach with respect 
of the linked yrast SD band in 152 Dy @, [H the spin of 
the initial state of this band has been suggested to be 

Io = 16. 5ft: this value has been used in the plotting 
of experimental kinematic moment of inertia in Fig. 5b. 
However, one cannot completely exclude also the value 
Io = 18.5ft (see discussion in Ref. dJ); this would in- 
crease experimental values bringing them closer to 
calculated ones. The rotational properties of the bands 
which come close to their maximum spin I max (as de- 
fined from the distribution of particles/holes over low- 
and high-j orbitals) are usually analyzed employing exci- 
tation energies plotted versus rigid rotor reference. This 
type of analysis is used for the bands in 109 Sb and 74 Kr, 
see Figs. 5c and 5d. One can see that at high spin above 
/ ~ 20ft, where the pairing is expected to be weak, the 
CRMF results describe the excitation energies very well. 
The band in 109 Sb is a classical example of a smoothly 
terminating band; it starts at low spin at near-prolate 
shape, then with increasing spin its shape becomes more 
and more triaxial and finally at I max = 41.5ft it termi- 
nates in a non-collective oblate state Q. This process 



of smooth band termination is associated with the grad- 
ual loss of collectivity which is nicely seen in the drop 
of the transition quadrupole moment Q t with increasing 
spin (see Fig. 5g). This feature is nicely reproduced in 
the calculations. On the contrary, the ground state band 
in 74 Kr reveals a new phenomenon of non-termination 
of rotational bands (see Ref. 58] for details): it remains 
collective at I ma x = 32ft. Contrary to the case of the 
band in 109 Sb, its transition quadrupole moment is more 
constant as a function of spin (see Fig. 5f). 

One can see that the CRMF calculations with the pa- 
rameter set NL3* provide a very successful description 
of different aspects of rotating nuclei. In general, the 
results of the calculations are very close to the ones ob- 
tained earlier with the NL3 parameterization of the RMF 
Lagrangian. 

As mentioned above, the density dependence of the 
parameter sets NL3 and NL3* is determined by the non- 
linear meson couplings in the isoscalar channel a. The 
vector fields oj and p, in particular, the isovector channel 
described by the p-meson are still coupled linearly and 
this leads to relatively large neutron radii and a very stiff 
equation of state for neutron matter [ijj as well as very 
stiff symmetry energies as a function of the density [4J] . 
Models with non-linear couplings in the p-channel as 
for instance the parameter set FSUGold introduced by 
Piekarewicz in Ref. [6~ij can improve this, however at 
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the cost of additional parameters. Density dependent 
covariant functionals eliminate the non-linear coupling 
terms and use, instead of that, at the cost of additional 
parameters, density dependent coupling constants 
g u {p), and g p {p) which depend on the baryon density p. 
This concept is evidently closer to the idea of density 
functional theory. On the other hand, one has to keep in 
mind that a proper treatment of the effects of vacuum po- 
larization by a renormalization procedure would require 
countcrterms which are nonlinear in the meson fields with 
the same polynomial structure [H, [l!| [2(| HH • In addi- 
tion, because of the relatively small surface diffuseness 
in nuclear systems and the narrow range, in which the 
density is changing from its value at saturation to zero, 
calculations of matrix elements with density dependent 
coupling constants and of the corresponding rearrange- 
ment terms require a considerably larger numerical effort, 
in particular, in the case of axially symmetric and triaxial 
applications. Considering all these points the non-linear 
version of RMF theory is still a very useful concept be- 
cause it combines simplicity with excellent predictions for 
many nuclear properties. 

In conclusion, in the present work we have reconsidered 
the well known parameter set NL3 after ten years and by 



a new fit with modern experimental data we introduced 
an improved parameterization NL3* for the RMF model, 
which contains only six phenomenological parameters. It 
is able to improve the description of nuclear masses and 
to cure some small problems observed previously with 
the NL3 force. At the same time, it provides excellent 
results for collective properties of vibrational and rota- 
tional character. 
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